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Abstract 

The  use  of  autoregressive  conditional  heteroskedastic  (ARCH)  models  to 
estimate  time-varying  hedge  ratios  suggests  that  conventional  procedures  may  not 
provide  appropriate  optimal  hedge  ratios.  However,  initial  results  using  ARCH 
models  raise  several  questions  regarding  time-varying  hedge  ratios.  First,  there 
are  several  alternative  representations  of  the  bivariate  generalized  ARCH 
(BGARCH)  model,  and  hedge  ratio  estimates  could  be  very  sensitive  to  model 
specifications.  Second,  the  available  models  should  be  tested  thoroughly  before 
selecting  an  appropriate  model.  Third,  alternative  approaches,  such  as  random 
coefficient  model,  to  capture  the  time  varying  nature  of  the  hedge  ratios  should 
be  explored.  Finally,  hedging  performance  in  terms  of  variance  reduction  of 
returns  with  alternative  procedures  is  an  important  practical  consideration. 
Using  data  on  the  cash  and  futures  price  changes  of  corn  and  soybeans,  following 
model  specification  tests,  a  diagonal  vech  parameterization  is  found  appropriate 
for  both  corn  and  soybeans.  The  diagonal  vech  parameterization  of  the  BGARCH 
model  provides  the  largest  reduction  in  the  variance  of  the  returns  portfolio 
which  is  consistent  with  diagnostic  and  specification  test  results.  For  both 
corn  and  soybeans,  the  constancy  hypothesis  of  the  hedge  ratio  was  rejected  when 
tested  against  the  random  coefficient  autoregressive  (RCMAR)  model.  However,  the 
RCMAR  procedure  does  not  lead  to  large  reductions  in  the  variance  of  returns. 
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1.     Introduction 

To  reduce  the  risks  of  price  uncertainty,  hedgers  often  use  futures 
contracts  to  cover  a  portion  of  their  spot  position.  As  the  economic  environment 
becomes  more  uncertain  and  volatile,  hedgers  need  more  precise  procedures  for 
identifying  appropriate  hedge  ratios.  Numerous  approaches  are  available  to 
estimate  hedge  ratios.  Traditionally,  ordinary  least  square  (OLS)  regression  of 
the  spot  price  on  the  futures  price  is  run,  with  the  slope  coefficient  being  the 
hedge  ratio  [e.g.  Ederington  (1979);  Anderson  and  Danthine  (1981);  Kahl  (1983); 
Benninga,  Eldor  and  Zilcha  (1984);  Carter  and  Loyns  (1985);  Brown  (1985)]. 
However,  this  procedure  is  statistically  inefficient  because  it  ignores  the 
heteroskedasticity  often  encountered  in  price  series  [e.g.  Kahl  (1984);  Bera, 
Bubnys,  and  Park  (1987);  Park  and  Bera  (1987);  Myers  and  Thompson  (1989)]. 

Myers  and  Thompson  (1989)  proposed  a  more  generalized  method  of  optimal 
hedge  ratio  estimation  which  incorporates  the  conditional  information  of  the 
mean,  variance  and  covariance  of  futures  and  spot  prices.  Two  special  problems 
still  require  attention.  One  is  the  constancy  of  hedge  ratio;  the  other  is  the 
regression  tendency.  When  either  the  variances  or  the  covariances  of  futures  and 
spot  prices  is  time-variant,  a  constant  hedge  ratio  is  inappropriate.  Also,  in 
the  context  of  a  market  model,  Blume  (1975)  observed  that  systematic  risk  (the 
beta  coefficient)  changed  over  time  in  a  regular  pattern.  Blume' s  empirical 
findings  identify  a  regression  tendency  of  the  beta  coefficient,  i.e.,  the 
estimated  beta  appears  to  take  less  extreme  values  and  exhibits  a  tendency  toward 
the  mean  [see  also  Bera  and  Kannan  (1986)]  .  We  might  expect  a  similar  phenomenon 
for  hedge  ratios,  when  its  OLS  estimate  obtained  from  historical  data  is  used  as 
the  ex  ante   hedge  ratio.   To  see  this,  we  plot  the  OLS  estimates  of  hedge  ratios 


for  corn  and  soybeans  in  Figures  1  and  2.  These  moving  estimates  are  obtained 
using  40  observations,  revising  the  data  set  through  time  by  adding  a  new 
observation  and  deleting  the  oldest  observation  each  period.^  The  figures 
clearly  display  regression  tendencies  over  time.  Thus,  the  constant  hedge  ratio 
generated  from  the  OLS  method  may  be  a  poor  estimate  for  one-time-period-ahead 
prediction. 

Recently,  with  the  development  of  autoregressive  conditional  hetero- 
skedastic  (ARCH)  [Engle  (1982)  and  generalized  ARCH  (GARCH)  Bollerslev  (1986)] 
models,  several  studies  estimating  time-varying  hedge  ratios  (THR)  have  been 
conducted.  The  time-varying  joint  distribution  of  cash  and  futures  price  changes 
has  been  examined  for  hedging  financial  instruments  (Cecchetti,  Cumby,  and 
Figlewski,  1988).  Bivariate  GARCH  (BGARCH)  models  also  have  been  applied  to 
estimate  the  THR  in  commodity  futures  [Myers  (1988);  Baillie  and  Myers  (1991); 
Myers  (1991)]  and  in  foreign  exchange  futures  [Kroner  and  Sultan  (1990)].  These 
recent  studies  suggest  that  conventional  hedging  procedures  are  inappropriate  and 
lead  to  hedge  ratios  which  differ  considerably  from  those  obtained  from  the  time- 
varying  framework.  The  initial  results  obtained  from  using  GARCH  models,  however 
raise  several  issues.  First,  there  are  several  alternative  specifications  of  the 
BGARCH  model,  and  hedge  ratio  estimates  could  be  very  sensitive  to  model 
specifications.  Second,  alternative  models  should  be  tested  thoroughly  before 
selecting  an  appropriate  model.  Third,  we  should  explore  other  alternative 
approaches,  such  as  random  coefficient  models,  which  may  capture  the  time-varying 
natures  of  hedge  ratios.  Finally,  hedging  performance  in  terms  of  variance 
reduction  of  returns  using  different  procedures  is  an  important  consideration. 
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The  data  set  will  be  discussed  in  Section  3 


In  this  paper,  we  address  the  above  issues  and  in  particular,  we  focus  on 
model  specification  and  empirical  comparison  of  BGARCH  and  random  coefficient 
models  for  the  hedge  ratio  estimation.  In  next  section,  the  hedging  rule  is 
derived  using  the  mean-variance  model.  In  Section  3,  different  BGARCH  models  and 
corresponding  estimation  results  are  presented.  We  also  conduct  diagnostic  tests 
to  identify  the  appropriate  model  specification.  In  Section  4,  we  apply  a  random 
coefficient  model  for  time -varying  hedge  ratios  and  examine  the  constancy  of 
hedge  ratio  represented  by  the  slope  coefficient.  We  compare  the  hedging 
effectiveness  of  the  random  coefficient  and  ARCH- type  models,  with  traditional 
OLS  hedging  approaches,  and  the  results  are  reported  in  Section  5.  Finally, 
Section  6  offers  some  concluding  remarks.  Our  analysis  is  similar  to  that  of 
Baillie  and  Myers  (1991);  however,  there  are  some  distinct  differences.  The  data 
sets  are  different,  and  we  use  nearby  futures  contract.  Our  use  of  random 
coefficient  model  in  the  context  of  hedge  ratio  estimation  is  new.  We  consider 
three  different  versions  of  the  BGARCH  model.  Lastly,  all  our  models  have  been 
subjected  to  various  diagnostic  checks,  specification  tests  and  model  selection 
procedures . 

2.     Time -varying  hedging  rule 

Define  f^  and  s^  as  the  changes  of  the  futures  and  the  spot  prices  between 
time  t  and  t-l,  respectively.  The  return  R^  on  the  investor's  portfolio  of 
holding  one  unit  of  the  spot  commodity  and  b^-i  units  of  the  future  contract  at 
time  t-l    can  be  written  as 

Rt  =  St  +  bfi-  ff  (1) 


In  the  absence  of  transaction  costs,  s^,  is  the  return  from  holding  the  cash 
position  from  c-l  to  t,  and  similarly  f^  is  the  return  from  holding  the  futures 
position  between  t  and  t-1.      We  rewrite  equation  (1)  in  terms  of  returns  as 

R,  =  R\  +  b,.i.  R^.  (2) 

Assume  that  the  investor  maximizes  the  utility  from  his  portfolio  by  choosing 
optimal  one-period  holdings  of  the  futures  contract  at  each  time  t.  The 
investor's  utility  function  is  assumed  to  be  monotonic,  continuous,  and  concave. 
Within  the  mean-variance  model,  the  expected  return  conditional  on  the  available 
information  at  C-l    is 

E   (i?Jij/,.i)  =  E   iRt"\^t-i^    +  V,£'(i?/|ijf,_,)  (3) 

where  Vt-i  i-S  the  information  set  at  time  C-l.  The  variance  of  returns  is  given 
by 

l^ar(RjVt-i)  =  ^^r   (RMVt-i)  +  bVi^ar(R^t  I  V't-i) 

+   2\.^Cov(R\,    R^|Vt-i)  (4) 

In  this  case,  solving  the  mean-variance  model  is  equivalent  to  maximizing  the 
expected  utility  function  which  in  turn  is  equivalent  to  maximizing 

E  (RMVt-i)  +  bt-iE(R^JVt-i  )  -  rt(I/ar(RMVt-i)  +  h\.^Var   (R^tl^t-i) 
-f-  2bt-i  Cov   (R\,    R^tlVt-i))  (5) 

where  r^.  is  the  risk  aversion  parameter.  To  obtain  the  utility  maximizing  hedge 
ratio,  differentiate  (5)  with  respect  to  bfi,  which  gives  the  optimal  number  of 
futures  contract  for  one  unit  of  spot  holdings.  The  utility  maximizing  hedge 
ratio,  b*t-i|V't-i»  is  given  by 

E(Rfj0t-i)  -  2r^cov    (R\,    R^Vt-i) 

bVilV-t-i  =  (6) 

2r^Var   (RMV-t-i) 


^ 


Assuming  E(ft|V't-i)  "  ^t-i.  that  is,  the  futures  market  is  unbiased,  E(R^tlV't-i)  is 
zero.   Then,  the  one-period  ahead  hedge  ratio  at  time  C-1    simplifies  to 

-Cov(/?/,i?/|Y,.,) 
jt>c-i  ~   :; (  '  ) 

This  time-varying  hedge  ratio  is  similar  to  the  conventional  hedge  ratio  except 
that  conditional  variance  and  covariance  replace  the  unconditional  variance  and 
covariance.  Because  conditional  moments  can  change  as  the  information  set  is 
updated,  the  hedge  ratios  also  can  change  over  time. 


3.     Time-varying  Hedge  Ratios  Estimation  Using  BGARCH  Models 


3.1    Specification  of  Selected  BGARCH  Models 

The  motivation  behind  using  BGARCH  models  in  the  context  of  hedge  ratio 
estimation  is  that  daily  commodity  future  and  spot  prices  react  to  the  same 
information,  and  hence,  have  non-zero  covariances  conditional  on  the  available 
information  set.   We  specify  the  model  as 


Re    =  V^s  ^  ^st  (8) 


e,   I  T,_i  ~  BN  {0,H^)  (9) 

where  R^t.  R^t  riow  denote  differences  of  logarithm  of  spot  and  futures  prices  from 
time  C    to  c-1    respectively,  6^=  (^st.  ^ft)  '  >   *t-i  is  the  information  set  up  to 


time  C-1,    BN  denotes  bivariate  normal  distribution,  and  H^  is  a  time-varying  2x2 
positive  definite  conditional  covariance  matrix.  Two  model  specifications  in  (8) 
are  used.''   Model  I  sets  /is=/if-0 ,  i.e.,  there  is  no  drift  term,  while  Model  II 
estimate  drift  parameters.   A  somewhat  general  form  of  H^  for  a  BGARCH  (p,q) 
model  can  be  written  as, 

(J  p 

vechiH^)    =  vechiC)    +  Y)  r_^vec^( 6^.16^-1)  +  Y  D^vechiH^.^)  (10) 

where  C  is  a  2x2  positive  definite  symmetric  matrix  and  F^  and  D^  are  3x3 
matrices.  The  "vech"  operation  stacks  the  lower  triangular  elements  of  a 
symmetric  matrix  in  a  column. 

The  parameterization  given  in  (10)  is  difficult  to  estimate  since  positive 
definiteness  of  H^  is  not  assured  without  imposing  nonlinear  parametric 
restrictions.  Moreover,  the  model  contains  too  many  parameters,  e.g. ,  for  p=q=l, 
Ht  has  21  parameters.  Here,  we  examine  several  specifications  of  U^.  One  simple 
assumption  that  could  be  made  is  to  specify  that  a  conditional  variance  depends 
only  on  its  own  lagged  squared  residuals  and  lag  values.  The  assumption  amounts 
to  taking  F  and  D  matrices  to  be  diagonal.  In  that  case,  vech(Ht)  of  a 
BGARCH(1,1)  model  is  given  by 


^  Other  specifications,  such  as  including  a  dummy  variable  for  daily  effect 
and  seasonal  dummies,  were  tried.  However,  estimates  for  dummy  variables  were 
not  significant. 


vech{H^)    = 


^ss,t 

Cs 

Jss      0 

0 

^sf,t 

= 

-sf 

+ 

0      7s/ 

0 

^ff,t 

0 

0 

h 

0         0 

2 

ss, t-1 

yff_ 

0     s^^ 

0 

' 

0 

0 

^ff_ 

h 

2 

ff,t-l 

'^ss,  t-1 

■,t-l^f,t-l 


ff,t-l 


(11) 


This  form  is  called  the  "diagonal  vech"  representation  of  H^^.   The  necessary 
conditions  for  this  H^  to  be  positive  definite  are 


c_j  >  0,  Cf  >  0,    c^Cf  -  Csf  >  0, 


(12) 


and  Y^5  >  0,  Y/f  >  0,  Y^^Y/f  "  Ysf  >  0. 
Baba,  Engle ,  Kraft,  and  Kroner  (1990)  suggested  another  parameterization  which 

is  almost  guaranteed  to  be  positive  definite.    This  form,   the  "positive 


definite,"  parameterization  is  written  as 


Ht  = 


'^ss    '^sf 
^fs    ^ff 


^ss    ^sf 
^fs    ^ff 


V 


yss  ysf 

7fs  7ff 


H^- 


t-1 


^ss    ^sf 
^fs     ^ff 


^ss,t-l  ^f,t-l^f,t-l 

2 

f,t-l^s,t-l  ^ff,t-l 


7ss  7s/ 
7fs  7/f 


(13) 


Note  that  the  number  of  parameters  to  be  estimated  for  this  specification  is  11. 

When  the  parameters  in  the  variance  covariance  function  are  time  invariant,  H^ 

becomes  a  constant  conditional  covariance  proposed  by  Myers  and  Thompson  (1989) , 

Bollerslev  (1990)  introduced  another  attractive  way  to  simplify  Hf   He  assumed 

that  the  conditional  correlation  between  £5^  ^^"^  ^ft  i^  constant  over  time  and 
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^t= 


33 


^3f 


^fS     ^ff 


(14) 


expressed  H^  as 


Ht    = 


,2     ,2 
^fs,t   ^ff,t 


's,  t 
0  h 


0 


1  Psf 
Psf      1 


's,  t 

0   h 


0 


(15) 


where  p^^  is  the  time- invariant  correlation  coefficient,  and  the  individual 
variances  hf  <.  and  h|  ^  ^^^^  assumed  to  be  standard  univariate  GARCH  process,  for 
example , 


It  is  clear  that  the  "constant  correlation"  representation  involves  only  7 
parameters.  Also,  positive  definiteness  of  the  specification  is  assured  if 
hg  t>0  hf  t>0  •  Because  of  these  attractive  features,  most  of  the  applications  of 
the  BGARCH  model  use  this  representation  [see  for  example,  Bollerslev  (1990), 
Baillie  and  Bollerslev  (1990),  and  Kroner  and  Sultan  (1990)]  .  However,  constancy 
of  correlation  is  a  very  strong  assumption  and  validity  of  model  (15)  remains  an 


empirical  question. 


3.2   Data  and  Selected  Sample  Statistics 

Daily  data  of  cash  and  futures  prices  for  corn  and  soybeans  are  used  for 
the  period  from  October  1988  to  December  1989.  The  cash  prices  are  central 
Illinois  elevator  bids.  March,  July,  November  nearby  futures  contracts  for 
soybeans  and  March,  July,  December  nearby  futures  contracts  for  corn  are  used. 
The  time  period  of  one  year  used  in  this  study  starts  from  the  beginning  day  of 
specific  futures  contract,  and  switches  to  the  nearby  contract  on  the  first  day 
of  expiration  month  for  the  specific  contracts.  The  number  of  observations  is 
311.  The  return  on  the  futures  and  cash  markets  are  defined  as  a  percentage 
change  of  each  prices,  i.e, 

R'\   =  1000  log  (S\/SVi) 
R'\   =  1000  log  (F\/F\.i) 
where  S^,^  and  F^j.  refers  to  i^*^  cash  and  futures  prices  respectively,  i  is  either 
corn  or  soybeans. 

Several  preliminary  diagnostic  tests  on  the  return  series  of  the  data  are 
conducted.  Using  the  Phillip  and  Perron  unit  root  test,  nonstationarity  of  the 
return  series  is  examined.  As  shown  in  Table  1,  the  unit  root  hypotheses  for  the 
all  return  series  of  cash  and  futures  of  corn  and  soybeans  are  rejected.  The 
sample  characteristics  of  univariate  unconditional  distribution  of  the  return 
series  of  corn  and  soybeans  are  summarized  in  Table  2.  Not  surprisingly,  cash 
and  futures  returns  for  the  corn  and  soybeans  show  typical  fat- tail  non-normal 
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distributions.  In  Table  3,  test  results  on  autocorrelation  of  the  series  and 
bivariate  normality  are  reported.  The  Ljung  and  Box  Q  statistic  is  calculated 
to  examine  the  autocorrelation  of  univariate  return  series.  There  is  no  evidence 
of  serial  correlation.  However,  The  Q^  statistics,  based  on  squared  return 
series,  provide  strong  indication  of  conditional  heteroskedasticity . 

Bivariate  normality  is  examined  using  the  Bera-John  (1983)  omnibus  test 
statistic 


C,   =  Ni'^  t!/S    +  T  ^'^ii   -   3)V24),  (16) 

i-l  i-l 

where  N  is  sample  size  and  Ti  =  ItiYilt/N,  Tii  =  It[yi]t/N,  where  y  =R  (X  -  m)  , 
X  is  the  vector  of  individual  observation  of  a  given  return  series  ,  m  is  the  mean 
of  X,  R  is  the  matrix  whose  eigen  vectors  are  the  same  as  the  eigen  vectors,  but 
whose  eigen  values  are  the  eigen  values  raised  to  -1/2,  of  the  sample  covariance 
matrix  of  X's.  Note  that  in  the  univariate  case,  T^  is  the  measure  of  skewness, 
and  Tii  is  the  measure  of  kurtosis.  Under  the  null  hypothesis  of  bivariate 
normality,  C3  is  asymptotically  distributed  as  x^  with  4  degrees  of  freedom. 
Given  the  high  values  of  the  test  statistics  bivariate  normality  is  rejected 
soundly  for  both  corn  and  soybeans  return  series. 
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3.3   Estimation  Results 

All  parameters  in  various  BGARCH  models  are  estimated  using  the  maximum 
likelihood  method.   The  log  likelihood  function  of  BGARCH  model  is  given  by 


i(e)  =  ^  Eic(0) 


C«l 


where      1,(6)    =  -^log  (27:) --|log|//,|-^e't^;'e(. 


(17) 


where  9  is  the  underlying  parameter  vector.  The  form  of  H^  depends  on  which 
specification  of  the  conditional  variance  is  used.  As  discussed  earlier,  three 
specifications  of  H^.,  namely  positive  definite,  diagonal  vech,  and  constant 
correlation  forms  of  BGARCH  are  estimated.  As  in  most  of  the  applied  work  in 
ARCH,  we  use  the  Berndt,  Hall,  Hall  and  Hausman  (1974)  algorithm  (BHHH)  for 


maximizing  1(S)  .        Starting  from  estimates  of  the  r   iteration,  the  (r+1) 
iteration  step  of  the  BHHH  algorithm  can  be  written  as 


:th 


Q(rn)  =    Q(r)  + 


f(^)M^ 


er  ae 


(18) 


where  the  first  derivatives  are  evaluated  at 


(r) 


Table  4  presents  the  estimation  results  of  the  constant  conditional 
variance  model.  The  estimation  results  of  BGARCH  models  with  various  specifi- 
cations are  reported  in  Tables  5,  6  and  7.  Using  likelihood  ratio  (LR) 
principle,  the  LR  statistics  are  presented  in  Table  8.  Under  the  null  hypothesis 

of  conditional  homoskedasticity ,  the  LR  statistics  are  asymptotically  distributed 

12 


as  x^(p) .  where  p  is  the  number  of  extra  parameters  in  different  BGARCH  models. 
From  the  results  of  Table  8,  it  is  clear  that  conditional  homoskedasticity  is 
rejected  overwhelmingly  when  tested  against  any  of  the  conditional  hetero- 
skedastic  models. 

Now,  considering  the  estimates  of  the  three  heteroskedastic  models,  we 
first  note  that  in  all  models  the  drift  terms  are  not  significant.  Also,  using 
the  LR  test,  the  presence  of  drift  terms  in  all  models  are  tested  under  the  null 
hypotheses  of  iJ,^=iJ,f=0  in  model  (II)  ,  which  has  a  asymptotic  x^(2)  distribution. 
The  null  hypothesis  are  not  rejected  based  on  LR  test  statistics,  0.04,  3.56  and 
3.1  for  corn  and  0.48,  1.72  and  5.15  for  soybeans,  respectively,  in  the  three 
BGARCH  model  (x^(2)=5.991  at  5%  significance  level).  Most  of  the  variance 
parameters  of  constant  correlation  and  diagonal  vech  models  are  significant  and 
values  of  the  estimated  parameters  seem  to  be  quite  reasonable.  For  the  diagonal 
vech  model,  the  conditions  for  positive  definiteness  of  H^  as  stated  in  (12)  are 
satisfied.  Also,  we  note  that  for  this  model,  the  estimates  of  y^^  +  6^^  and  y^f- 
+  6ff  are  all  close  to,  but  less  than,  one.  On  the  other  hand,  in  the  positive 
definite  parameterization  model,  many  of  the  parameter  estimates  are  not 
significant.  Also,  for  this  model,  the  sum  of  the  coefficients  of  lagged  squared 
residuals  and  past  conditional  variance  sometimes  exceeds  one.  For  example,  in 
the  case  of  corn,  7f£  +  Sff  takes  on  values  of  1.067  and  1.071,  respectively, 
for  Models  I  and  II.   This  might  lead  to  instability  of  unconditional  higher 
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moments.  Of  course,  it  is  not  possible  to  choose  a  model  just  on  the  basis  of 
significance  and  implications  of  parameter  estimates.  For  model  selection,  we 
subject  three  BGARCH  models  to  more  rigorous  model  testing.  In  the  next  section, 
we  perform  a  number  of  diagnostic  checks,  such  as,  conditional  bivariate 
normality,  auto-correlation  and  conditional  heteroskedasticity  tests,  using  the 
standardized  residuals  obtained  from  the  estimated  models.  Model  selection 
criteria  such  as  Akiake's  Information  Criteria  (AIC)  and  Baysian  Information 
Criteria  (BIC)  are  also  used. 


3.4   Model  Specification  Tests 

Economic  theory  does  not  provide  specific  guidelines  for  the  appropriate 
parameterization  of  the  BGARCH  model.  Therefore,  the  selection  of  an  appropriate 
model  is  essentially  an  empirical  question.  Based  on  the  diagonal  vech  and 
positive  definite  parameterizations ,  we  conduct  the  Ljung-Box  Q,  Q^,  and 
conditional  normality  tests  of  the  standardized  residuals,  H^'^''^  e^.  For  the 
constant  correlation  BGARCH  model,  the  results  are  reported  in  Table  9.  As  in 
the  original  data  from  the  Q  statistic  values,  we  note  that  there  is  no  evidence 
of  serial  correlation  in  the  residuals.  The  Q^  statistic,  based  on  squared 
residuals,  provides  an  indirect  test  of  conditional  heteroskedasticity.  None  of 
the  Q^  statistics  are  significant,  which  contrasts  sharply  with  the  Q^  values  for 
the  return  series  reported  in  Table  3.    This  indicates  that  the  constant 

14 


correlation  model  does  take  account  of  the  conditional  heteroskedasticity  present 
in  the  data.  Finally,  the  Bera-John  test  statistics  for  bivariate  normality  are 
much  lower  than  those  for  the  original  series;  for  instance,  for  the  original 
soybean  series,  the  test  statistic  was  1047.75.  This  was  reduced  to  206.27  for 
the  standardized  residuals  using  Model  II.  However,  all  the  statistics  are  still 
significant  indicating  that  constant  correlation  BGARCH  does  not  take  into 
account  all  the  nonnormality  present  in  the  data.  Again,  the  main  source  of 
nonnormality  is  the  high  kurtosis  values. 

The  diagnostic  test  results  for  the  diagonal  vech  and  positive  definite 
BGARCH  models  are  reported  in  Tables  10  and  11,  respectively.   All  the  skewness 
and  kurtosis  coefficients  indicate  that  the  conditional  distributions  are  still 
fat-tailed  non-normal.     The  calculated  Bera-John  test,  rejects  the  null 
hypothesis  of  bivariate  conditional  normality  for  both  the  corn  and  soybean  data.- 
However,  the  use  of  BGARCH  filter  results  in  a  significant  reduction  in  the  non- 
normality.   For  example,   the  value  of  the  Bera-John  statistic,  563.16,  for  the 
unconditional  distribution  of  return  series  of  corn  is  reduced  to  87.61  and 
34.44,  respectively,  for  the  diagonal  vech  and  positive  definite  representations 
with  no  drift  terms.   The  Q(12)  statistics  indicate  that  all  the  specifications 
have  no  autocorrelation  except  model  (II)  under  the  diagonal  vech  parameter- 
ization for  corn.   However,  except  for  a  few  cases,  the  Q^(12)  statistics  are 
significant. 

15 


0 


In  order  to  select  an  appropriate  specification  of  BGARCH  model,  AIC  and 
BIG  also  are  calculated,  and  are  given  in  Tables  5,  6  and  7.  Model  I  for  corn 
and  soybeans  in  the  positive  definite  parameterization  have  the  smallest  AIC  and 
BIG  among  possible  specifications.  Based  on  specification  tests,  model  selection 
criteria  and  our  earlier  critical  comments  on  the  parameter  estimates  of  the 
positive  definite  form,  the  diagonal  vech  model  seems  to  be  appropriate  for  the 
time-varying  hedge  ratios  estimation  for  both  corn  and  soybeans. 

The  BGARGH  models  provide  the  time-varying  conditional  variances  and 
covariances  of  R^^  ^^*^  ^^t-  Thus,  the  time-varying  hedge  ratio  at  time  t-1  can 
be  obtained  from, 

^\f  t 

bVi  =  '■ =  HRt  (19) 

The  time-varying  hedge  ratios  of  corn  and  soybeans  calculated  from  the 
above  equation  are  plotted  in  Figures  3,  4,  and  5,  respectively.  In  Figures  3 
and  4,  the  hedge  ratios  based  on  diagonal  vech  and  constant  correlation  BGARGH, 
and  constant  hedge  ratio  estimates  of  0.93  based  on  OLS  method  are  illustrated 
for  corn.  The  means  and  standard  deviations  of  the  time-varying  hedge  ratios  of 
constant  correlation  and  diagonal  vech  BGARGH  strategies  for  corn  are  0.970  and 
0.931,  and  0.101  and  0.125  over  the  sample  period.  The  estimated  time-varying 
hedge  ratios  from  the  two  specifications  have  similar  movements  around  its  mean. 
These  figures  are  quite  different  from  the  moving  OLS  estimates  of  the  hedge 
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ratio  in  Figure  1.  This  implies  attempts  to  obtain  time-varying  hedge  ratios 
using  OLS  with  most  recent  observations  may  not  be  very  fruitful.  In  the  case 
of  soybeans  (Figure  5),  the  constant  hedge  ratio  from  OLS  is  0.897.  The  mean  and 
standard  deviation  of  time-varying  hedge  ratios  of  diagonal  vech  BGARCH  strategy 
are  0.962  and  0.199  over  sample  period.  Unlike  the  case  of  corn,  for  soybeans, 
this  sample  mean  is  quite  different  from  the  OLS  hedge  ratio.  In  Figure  5,  the 
THR  has  one  substantial  drop  at  the  beginning  of  the  period  (October  -  November 
1988)  and  high  peak  around  July  1989.  These  large  changes  in  the  hedge  ratios 
occur  near  maturity  of  the  harvest  contracts  (October  -  November)  and  during  the 
summer  when  the  crops  are  most  susceptible  to  the  effects  of  weather. 


4.     Comparison  with  Random  Coefficient  Model 


The  conventional  hedge  ratio  model  is  given  by 

R,'  =  a   ^   Pi?/  +  €,  ^20) 

€j.  ~  i.i.d.  N{0,  a^) 

The  primary  objection  to  this  model  is  the  time  invariance  of  the  coefficient  /9 . 
In  the  earlier  section,  we  discussed  the  use  of  ARCH  models  to  capture  the 
dynamic  nature  of  hedge  ratios.  The  random  coefficient  model  (ROM)  is  another 
approach  to  take  account  of  the  time -varying  nature  of  hedge  ratios  and  in  this 
section  we  explore  the  possibility  of  using  RCM  as  an  alternative  to  ARCH  models. 
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In  a  different  context,  this  type  of  model  has  been  used  before  to  estimate  time- 
varying  coefficients.  For  example,  Bos  and  Newbold  (1984)  used  RCM  to  estimate 
time-varying  systematic  risk  in  the  market  model  framework. 

A  time-varying  coefficient  model  with  an  autoregressive  specification 
(denoted  RCMAR)  is  written  as 

^/  =  «  ^  Pt  ^/  -^  V, 

(P,  -]r)  =  (j)(p,.,  -]J)  +  fi, 

where  ^^  is  a  time-varying  coefficient,  and  i/^  and  /i^  are  independent  and 
identically  distributed  (IID)  random  variables  with  means  zero  and  variances  a^^ 
and  a^^  respectively,  and  |0|(l.  When  0  =  0,  the  RCMAR  is  identical  to  the 
standard  RCM  of  Hildreth  and  Houck  (1968).  The  parameters  of  the  model  can  be 
estimated  using  the  nonlinear  maximum  likelihood  procedures  described  in  Pagan 
(1980)  .  The  estimated  results  are  shown  in  Tables  12  and  13  for  corn  and 
soybeans,  respectively.  All  parameter  estimates  for  corn  and  soybeans  are 
significant  except  the  constant  term  a. 

The  most  important  issue  here  is  to  test  the  constancy  of  beta  coefficient. 
For  simplicity,  first  consider  the  case  of  RCM,  i.e.,  (f)  =  0 .  The  relevant 
testing  problem  is  to  test  H^ :  a^  =  0  against  Hg :  cr^  >  0 .  From  the  estimates  of 
RCM  models,  we  note  that  a^  is  highly  significant  for  both  corn  and  soybeans. 
Also,  the  asymptotic  x^(l)  LM  statistics  for  testing  a^  =  0   are  16.599  and  3.977 
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for  corn  and  soybeans,  respectively.    Thus,  the  null  hypotheses  of  fixed 
parameter  is  rejected  at  both  the  1  and  5  percent  significance  levels. 

When  4>  in  the  time-varying  coefficient  model  (21)  is  not  equal  to  zero, 
then  testing  Hg:  cr^=0  is  more  complicated.  If  the  beta  coefficient  is  not 
constant  over  time,  under  (21),  ^ar{fi^)~a^/l-4>^.  Under  the  null,  W^:  a^=0 ,  the 
parameter  (f)  is  not  identified.  Therefore,  a  test  can  not  be  conducted  using  the 
conventional  large  sample  tests  such  as  the  LR,  LM,  or  Wald  tests  because  the 
nuisance  parameter,  (f> ,  is  identified  only  under  alternative  hypothesis.  This 
phenomenon  results  in  the  violation  of  regularity  conditions  for  the  standard 
testing  procedures.  For  example,  in  the  formulation  of  the  LM  test,  the 
information  matrix  derived  from  the  likelihood  function  with  the  restrictions  of 
a^=0  becomes  singular.  Here,  we  follow  the  Davies  (1977)  approach  for  testing 
a^=0 .  This  procedure  involves  computing  the  standard  LM  statistic  for  a  given 
value  of  <j> ,    say  LM(0)  .   Then  the  test  is  based  on  a  critical  region  of  the  form 

{  SUP    LM{^)      >  c  )  /oo\ 

(J)G(-1,1)  ^'^^'' 

where  c  is  a  suitably  chosen  constant.  Davies  (1987)  provided  an  approximation 
to  the  p-value  of  the  test  [for  a  description  of  the  computation  of  the  p-value, 
see  Bera  and  Higgins  (1992)].  The  Supremum  of  LM((^)  for  corn  and  soybeans  are 
11.254  and  3.091  at  0=-O.76  and  0.28,  respectively.  [It  is  interesting  to  note 
that  the  corresponding  maximum  likelihood  estimate  of  (f)   in  Tables  12  and  13  are  - 
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0.758  and  0.28].  The  p-values  of  the  Davies  test  are  0.00  and  0.00099 
respectively  for  corn  and  soybeans.  Therefore,  the  null  hypothesis  of  constant 
coefficient  model  is  rejected  for  both  corn  and  soybeans. 

The  test  of  a  RCM  against  the  alternative  of  RCMAR  is  carried  out  with  a 
LR  test  under  the  following  testable  hypotheses 

Ho-'      <t>  =  0 

H^   :     4)  '^  0  ^^^^ 

The  calculated  LR  statistics  are  15.699  and  4.169  for  corn  and  soybeans 
respectively.  Thus,  the  null  hypothesis,  RCM,  is  rejected,  although  for  soybeans 
the  rejection  is  only  at  5  %  significance  level. 

These  test  results  indicate  that  the  RCMAR  model  provides  a  good  represen- 
tation of  the  data  set.  In  order  to  trace  time-variant  hedge  ratios,  i.e.,  /3t, 
we  employ  a  fixed  interval  smoothing  process  using  the  conditional  means  and 
variances  of  ^^  over  time  period,  t=l,2,..,T  [see  Newbold  and  Bos  (1985,  p. 39)] 


=  E[^,\R^,...R^^;R,,.  ..R,^]  (24) 

=  P(t|t)  +^jp(t^-i|r)    -P(t  +  l|f:)] 


Pt  =  P(t|r) 


=  Pit\t)    -  A^[P{t^l\t)    -  P{t  +  l\T)]A^ 
where    A^.  =  P{t\t)  (^P{t+l\t)  ~'^ . 


(25) 


The  smoothing  process  begins  by  setting  t=T-l  in  the  (24)  and  (25).  The  maximum 

likelihood  estimates,  y9 ,  <p ,    a^,    and  a^ ,    presented  in  Tables  12  and  13,  are  used 

for  the  fixed  coefficient  in  the  estimation  of  stochastic  parameter  over  sample 
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period.  The  plots  of  p^,  time-varying  hedge  ratios  for  corn  and  soybeans,  are 
given  in  Figures  6  and  7  respectively.  For  corn,  there  is  some  similarity 
between  the  plots  of  THR  using  RCMAR  and  those  reported  in  Figures  3  and  4  using 
BGARCH  model;  although  in  Figure  6  we  observe  more  fluctuation.  For  soybeans, 
quite  interestingly,  the  behavior  THR  obtained  from  BGARCH  and  RCMAR,  as  can  be 
seen  from  Figures  5  and  7  are  very  close. 


5.     Hedging  Performance  of  BGARCH  and  RCMAR 


The  results  of  the  diagnostic  tests  for  the  BGARCH  and  RCMAR  models  led  us 
to  identify  alternative  specifications  for  estimating  the  time-varying  hedge 
ratios.  However,  the  hedging  performance  of  the  models  still  remains  a 
practical  question.  Here,  the  in-  and  out-of -sample  hedging  performance  of  the 
BGARCH  and  RCMAR  models  are  investigated.  Under  the  assumption  that  the  time- 
varying  conditional  variance  and  covariance  estimated  from  the  diagonal  vech 
BGARCH  model  is  the  appropriate  generating  process,  the  percentage  reduction  in 
the  conditional  variance  of  the  returns  portfolio  relative  to  the  no  hedging 
scenario  is  used  to  evaluate  hedging  performance.  In-sample  results  are  based 
on  the  311  observations  used  in  the  estimation  process,  while  the  out-of -sample 
hedging  performance  is  assessed  over  the  next  60  observations . 
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The  results  of  the  in-  and  out-of -sample  hedging  performance  are  reported 
in  Tables  14  and  15,  respectively.  For  the  corn,  using  the  BGARCH  models,  the 
l|  in-sample  hedge  performance  is  very  similar,  with  an  average  variance  reduction 
of  79  percent.  For  soybeans,  the  BGARCH  representations  also  produce  similar 
reductions  in  the  variance  of  returns,  except  for  the  positive  definitive 
parameterization.  In  the  context  of  model  I,  which  most  accurately  represented 
the  structure  of  the  data,  the  diagonal  vech  specification  provides  the  largest 
variance  reduction,  79.42  and  77.00  percent  reductions  for  corn  and  soybeans, 
respectively.  This  is  consistent  with  the  results  of  the  diagnostic  and 
specification  tests  used  to  select  the  appropriate  conditional  heteroskedastic 
model.  For  both  corn  and  soybeans,  the  RCMAR  performs  rather  poorly,  reducing 
the  corn  and  soybean  variances  by  73.70  and  73.23  percent,  respectively,  which 
is  inferior  even  to  the  constant  hedging  (OLS)  procedure. 

Out-of -sample  results  verify  the  attractiveness  of  the  diagonal  vech 
parameterization;  its  use  leads  to  the  largest  variance  reductions  for  both 
commodities  and  models.  The  other  BGARCH  models  and  the  RCMAR  perform  very 
poorly,  highlighting  the  importance  of  correct  specification  procedures  for  out- 
of-  sample  use  of  the  BGARCH  models. 

Overall,  the  results  indicate  the  usefulness  of  appropriately  specified 
BGARCH  models  for  establishing  hedge  ratios.  These  findings  are  consistent  with 
Baillie  and  Myers  (1991)  but  somewhat  different  from  Myers  (1991)  who  encountered 
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little  difference  in  hedging  performance  between  time -varying  and  other 
procedures.  Clearly,  the  disparate  findings  may  be  a  result  of  differences  in 
the  commodities  and  data,  as  well  as  the  specification  procedures  employed  here. 
The  poor  results  associated  with  the  RCMAR  framework  caution  against  the  use  of 
a  procedure  that  does  not  directly  specify  and  estimate  the  changing  nature  of 
variance  matrix  which  is  fundamental  to  the  measurement  of  time -varying  hedge 
ratios . 

While  the  BGARCH  model  shows  potential  for  improved  risk  management, 
further  study  is  needed  to  assess  the  costs  of  its  specification  and  implemen- 
tation relative  to  the  gains  in  variance  reduction.  The  implementation  of  the 
BGARCH  framework  can  require  frequent  and  costly  position  changes  in  the  futures 
market.  In  addition,  the  estimation  and  continual  updating  of  models  for 
practical  use  can  be  time-consuming  and  costly.  Assessment  of  the  hedging 
performance  of  these  models  in  framework  which  explicitly  incorporates  these 
costs  will  provide  a  more  complete  understanding  of  the  usefulness  of  time- 
varying  procedures  for  managing  price  risk. 


6.     Conclusion  and  Further  Research 


Various  BGARCH  and  RCMAR  models  were  applied  to  estimate  time-varying  hedge 
ratios.   This  remedies  the  constancy  and  regression  tendency  of  hedge  ratios 
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based  on  conventional  OLS  estimation.  After  rigorous  tests  of  model  specifi- 
cation, diagonal  vech  parameterization  is  found  to  be  appropriate  for  both  corn 
and  soybean.  The  BGARCH  hedge  ratios  resulted  in  better  performance  of  reducing 
the  variance  of  the  portfolio  return  for  corn  and  soybeans.  A  diagonal  vech 
parameterization  provides  largest  reduction  of  the  variance  of  the  portfolio 
return.  This  is  consistent  with  the  results  of  diagnostic  and  specification  test 
for  the  representation  of  conditional  second  moment.  For  both  corn  and  soybeans, 
the  constancy  hypothesis  of  hedge  ratio  from  OLS  was  rejected  against  random 
coefficient  autoregressive  (RCMAR)  model.  However,  the  RCMAR  strategy  shows  poor 
performance  for  the  reduction  of  variance  of  return. 

There  are  a  number  of  issues  that  require  further  attention.  In 
particular,  the  relationship  between  the  BGARCH  and  random  coefficient  models  is 
an  interesting  issue.  For  our  data  set,  the  performance  of  RCM  is  not  very 
encouraging.  However,  from  a  theoretical  point  of  view,  it  is  a  viable 
alternative  to  the  BGARCH  model  for  estimating  time-varying  hedge  ratios,  and 
might  work  better  for  other  data  sets.  The  fundamental  difference  between  these 
two  models  is  that  the  BGARCH  models  consider  the  joint  conditional  distribution 
of  the  spot  and  future  prices  while  RCM  specifies  a  conditional  mean  model  for 
the  spot  price  given  the  future  prices.  In  RCM  the  time  varying  nature  of  the 
model  is  incorporated  through  the  mean  equation  and  in  BGARCH  models,  this  is 
achieved  through  the  variance  specification.    From  the  behavior  of  the 
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standardized  residuals,  we  noted  that  the  BGARCH  models  with  conditional 
normality  assumption  does  take  account  of  high  kurtosis  to  some  extent.  However, 
a  part  of  non-normality  still  remains  unexplained.  Therefore,  use  of  bivariate 
conditional  t-distribution  would  be  worthwhile. 
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Table  1.     Phillips-Perron  Test  for  Unit  Roots  in  Return  Series 

Corn  Soybean 


Cash  Futures                  Cash       Futures 

Z(t;)             -18.30  -17.39                 -19.60      -15.92 

Z($i)             167.98  151.82                 244.14       91.85 

Z(t;)             -18.77  -17.94                 -19.37      -17.22 

Z($2)             131.70  129.86                 216.40      115.92 

Z($3)             111.64  101.16                 210.76       43.07 

Note: 

Critical          Z(t;)  Z($i)       Z(t;)        Z($2)        2(^3) 

Values 

1%               -3.43  6.43        -3.96       8.27        6.09 

5%               -2.86  4.59        -3.66       6.25        4.68 

The  tests  are  based  on  the  following  regressions, 

y^  =  fl  +  P(t-r/2)  +  ay,_^  +  u^ 

Z(t„.)  test  for  H^  :                  a*=  1 

Z(OJ  test  for  H^   :   n'  =  0  a*=  1 

Z(  tg)  test  for  H^  :                    a  =  1 

Z($2)  ^^^^  -^<^-^  ^o   '       (1=0  5=1   j3=0 

Z($3)  test  for  H^  I                  a  =  1   P  =  0 
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Table  2. 

Sample  ; 

Statistics  of  t 

:he  Return 

Series  —  Jan. 

1989-Dec. 

1990 

Variable 

Mean 

Std.  Dev. 

Skewness 

Kurtosis 

Minimum 

Maximum 

Corn 

Spot 
Futures 

-.55555 
-.73153 

13.234 
12.503 

-.50831 
-.24675 

4.4112 
4.8365 

-56.05 
-50.06 

37.18 
43.34 

Soybeans 

Spot 
Futures 


1.1376 
1.0438 


14.390 
13.414 


61165 
37488 


4.2226 
4.0415 


-61.70 
-46.80 


36.29 
43.68 


Futures  Contracts 


Corn 
Soybean 


May,  July,  December 
March,  July,  November 
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Table  3.     Test  Results  of  Autocorrelation  and  Blvariate  Normality  for  the 
Return  Series 


Corn 


Soybean 


SPOT 


FUTURES 


SPOT 


FUTURES 


Q(12) 
Q'(12) 


10.84 
35.04 


14.91 
52.41 


19.83 
39.30 


10.25 
35.45 


Bivariate 
Normality 

Ti 

Til 


•0.38  -0.23 

4.40  9.38 

563.16 


■0.41  0.83 

3.87  11.75 

1047.75 


•11 


Skewness 
Kurtosis 
Bera- John's  LM  statistic, 


Critical  values  of  x  (p)  are: 


1% 

5% 

p  =  4 

13.277 

9.488 

p  =  12 

26.217 

21.026 
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Table  4.     Estimation  Results  of  the  Bivariate  Models  with  Constant 
Conditional  Covariance  Matrix 


I^^ 


Corn 


II 


2) 


Soybean 


1) 


Il2) 


0.634 
(0.824) 


0.989 
(0.833) 


/^f 


0.790 
(0.765) 


1.085 
(0.762) 


176.257 
(10.985) 


183.170 
(12.386) 


206.418 
(13.691) 


195.557 
(13.308) 


'sf 


153.328 
(3.001) 


153.392 
(3.231) 


161.755 
(5.073) 


152.747 
(5.324) 


168.182 
(11.106) 


163.519 
(10.835) 


180.204 
(11.547) 


174.737 
(11.323) 


Log 
Likelihood 


-2233.55 


2231.98 


2322.59 


2321.63 


1) :    Refers  to  the  model  specification  without  drift  term 
2) :    Refers  to  the  model  specification  with  drift  term 
(  ) :   standard  error 
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Table  5 


Estimation  Results  of  the  Constant  Correlation  BGARCH  Model 


Corn 

Soybean 

Ii> 

Il2) 

I^^ 

Il2) 

Ms 

0.076 
(0.747) 

-0.410 
(0.698) 

A^f 

0.114 
(0.691) 

-0.182 
(0.683) 

Cs 

23.936 
(6.808) 

23.305 
(7.017) 

43.937 
(15.341) 

44.155 
(15.152) 

Cf 

14.245 
(3.323) 

13.540 
(3.200) 

31.477 
(11.664) 

31.472 
(11.644) 

7ss 

0.136 
(0.0436) 

0.133 
(0.044) 

0.200 
(0.037) 

0.201 
(0.038) 

7ff 

0.121 
(0.019) 

0.117 
(0.019) 

0.133 
(0.039) 

0.133 
(0.039) 

^SS 

0.738 
(0.064) 

0.387 
(0.066) 

0.605 
(0.090) 

0.604 
(0.089) 

5ff 

0.790 
(0.022) 

0.799 
(0.022) 

0.703 
(0.079) 

0.702 

(0.079) 

Psf 

Log 

Likelihood 
AIC 
BIC 

0.887 
(0.011) 

-2209.02 
4432.04 
4458.22 

0.887 

(0.011) 

-2209.00 
4436.04 
4469.70 

0.857 
(0.010) 

-2294.14 
4602.28 
4628.46 

0.857 

(0.011) 

-2293.90 
4605.80 
4639.46 

1): 
2): 

(  ) 
AIC 
BIC 


Refers  to  the  model  specification  without  drift  term 

Refers  to  the  model  specification  with  drift  term 

standard  error 

Akaike  Information  Criteria 

Bayesian  Information  Criteria 
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Table  6.     Estimation  Results  of  the  BGARCH  Models  with  Diagonal  Vech 
Parameterization 


Corn 

S( 

Dybean 

Ii> 

Il2> 

1^' 

Il2) 

/^s 

0.121 
(0.780) 

0.668 
(0.701) 

/^f 

0.454 
(0.730) 

0.938 
(0.668) 

Cs 

68.418 
(16.299) 

63.880 
(15.283) 

58.103 
(10.183) 

62.705 
(11.023) 

Csf 

47.854 
(12.531) 

44.148 
(11.394) 

45.811 
(7.836) 

48.229 
(8.100) 

Cf 

37.956 
(10.575) 

33.362 
(8.834) 

37.862 
(7.079) 

38.246 
(7.017) 

7ss 

0.197 
(0.053) 

0.189 
(0.051) 

0.345 
(0.043) 

0.350 
(0.041) 

7sf 

0.145 
(0.43) 

0.143 
(0.040) 

0.314 
(0.043) 

0.310 
(0.040) 

7ff 

0.132 
(0.039) 

0.128 
(0.037) 

0.304 
(0.048) 

0.291 
(0.044) 

5ss 

0.454 
(0.102) 

0.479 
(0.099) 

0.485 
(0.048) 

0.463 
(0.047) 

5sf 

0.561 
(0.095) 

0.584 
(0.089) 

0.543 
(0.040) 

0.532 
(0.381) 

<5ff 

Log 

Likelihood 
AIC 

BIG 

0.641 
(0.079) 

-2193.17 
4404.34 
4437.99 

0.673 
(0.068) 

-2191.39 
4404.78 
4445.92 

0.589 
(0.043) 

-2181.97 
4381.94 
4415.60 

0.593 
(0.041) 

-2182.83 
4387.66 
4428.80 

1) :    Refers  to  the  model  specification  without  drift  term 
2):    Refers  to  the  model  specification  with  drift  term 


(  ) 
AIC 

BIG 


standard  error 

Akaike's  Information  Criteria 

Bayesian  Information  Criteria 
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Table  7.     Estimation  Results  of  the  BGARCH  Models  with  Positive  Definite 
Parameterization 


Corn 

Soybean 

I^^ 

Il2> 

I^^ 

Il2) 

Ms 

0.379 
(0.762) 

0.703 
(0.813) 

/^f 

0.693 
(0.742) 

0.954 
(0.758) 

Css 

29.921 
(13.936) 

29.462 
(13.768) 

50.222 
(20.064) 

51.671 
(30.470) 

Csf 

12.091 
(6.705) 

11.938 
(6.593) 

48.176 
(4.778) 

48.170 
(4.983) 

Cff 

11.129 
(4.206) 

10.824 
(4.382) 

47.140 
(25.029) 

46.028 
(26.575) 

7ss 

0.770 
(0.125) 

0.780 
(0.132) 

0.741 
(0.119) 

0.716 
(0.117) 

7sf 

0.160 
(0.076) 

0.158 
(0.079) 

-0.0008 
(0.119) 

-0.028 
(0.116) 

7fs 

-1.027 
(0.151) 

-1.046 
(0.161) 

-0.591 
(0.128) 

-0.547 
(0.126) 

7ff 

-0.479 
(0.101) 

-0.474 
(0.106) 

0.239 
(0.128) 

0.275 
(0.127) 

^SS 

0.341 
(0.212) 

0.347 
(0.212) 

0.595 
(0.933) 

0.620 
(0.093) 

^sf 

-0.012 
(0.165) 

-0.013 
(0.158) 

-0.142 
(0.099) 

-0.112 
(0.098) 

5fs 

0.507 
(0.206) 

0.499 
(0.206) 

0.239 
(0.079) 

0.207 
(0.084) 

5ff 

Log 
Likelihood 

AIC 
BIG 

0.918 
(0.153) 

-2175.02 
4372.04 
4413.17 

0.920 
(0.146) 

-2173.92 
4373.84 
4422.46 

0.951 
(0.686) 

-2166.10 
4354.20 
4395.34 

0.925 
(0.074) 

-2163.75 
4353.50 
4402.12 

1): 
2): 

(  ) 
AIC 
BIG 


Refers  to  the  model  specification  without  drift  term 

Refers  to  the  model  specification  with  drift  term 

standard  error 

Akaike's  Information  Criteria 

Bayesian  Information  Criteria 
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Table  8.     LR  Statistics  for  Testing  Constant  Covariance  Model  Against 
Different  BGARCH  Models 


Model 
I 

II 


Model 
I 

II 


Corn 


Constant  Correlation    Diagonal  Vech 


49.08 
(4) 

45.96 
(4) 


Constant  Correlation 

56.90 
(4) 

55.46 
(4) 


80.76 
(6) 

81.18 
(6) 

Soybeans 

Diagonal  Vech 

281.24 
(6) 

277.60 
(6) 


Postive  Definite 

117.06 
(8) 

116.12 
(8) 


Postive  Definite 

312.98 
(8) 

315.76 
(8) 


1)  Numbers  in  the  parentheses  are  the  degrees  of  freedom  for  the  LR  test, 

2)  Critical  values  of  x^(p)  ^^^  given: 


p  =  4 
p  =  6 
p  =  8 


1% 
13.277 
16.811 
20.090 


5% 
9.488 
12.592 
15.507 
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Table  9     Test  Results  of  Autocorrelation  and  Conditional  Normality  for  the 
Constant  Correlation  BGaRCH  Model 


Ii> 


Corn 


llZ) 


Soybean 


Ii> 


II 


2) 


Q(12) 


St 

'ft 


9.15 
14.19 


9.12 
14.17 


17.92 
10.72 


18.00 
10.80 


Q^(12) 


^^t 

11.57 

11.41 

^\t 

11.14 

11.05 

Conditional 

Normality 

Ti 

0.17 

0.14 

T2 

-0.63 

-0.64 

Til 

6.23 

6.25 

T22 

5.82 

5.87 

c. 

261.16 

265.70 

6.59 
11.30 


6.75 
11.39 


-0.06 

-0.12 

-0.60 

-0.61 

6.99 

6.51 

4.21 

4.45 

243.21 

206.27 

1) :    Refers  to  the  model  specification  without  drift  term 
2) :    Refers  to  the  model  specification  with  drift  term 
Ti :    Skewness,  i=l ,  2.  l=spot  returns,  2=futures  returns, 
Ti^:   Kurtosis,  i=l ,  2.  l=spot  returns,  2=futures  returns, 
C3:    Bera-John's  LM  statistic,  x^i^) ■ 
Critical  values  of  x^(p)  are: 


1% 

5% 

p   =      4 

13.277 

9.488 

p  =   12 

26.217 

21.026 
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Table  10.    Test  Results  of  Autocorrelation  and  Bivariate  Conditional  Normality 
in  the  BGARCH  Models  with  Diagonal  Vech  Parameterization 


Corn 


I^) 


II 


2) 


Soybean 


I^^ 


II 


2) 


Q(12) 


fist 

€ft 


19.08 
11.28 


24.24 
17.86 


12.31 
8.44 


12.58 
6.05 


Q^(12) 


«^t 

7.19 

44.35 

^\t 

19.44 

3.41 

Conditional 

Bivariate 

Normality 

Ti 

0.054 

0.33 

T2 

-0.24 

-0.55 

Tn 

5.54 

6.44 

T22 

3.19 

4.15 

c. 

87.61 

191.89 

12.37 
29.73 


12.30 
23.03 


0.057 

-0.032 

-0.45 

-0.47 

5.11 

5.09 

3.67 

3.59 

74.37 

72.84 

1) :  Refers  to  the  model  specification  without  drift  term 

2) :  Refers  to  the  model  specification  with  drift  term 

Ti :  Skewness,  i=l ,  2.  l=spot  returns,  2=futures  returns, 

Tii!  Kurtosis,  i=l ,  2.  l=spot  returns,  2=futures  returns, 

C3:  Bera-John's  LM  statistic,  x'^(^) - 


Critical  values  of  x  (p)  are 


1% 

5% 

p  =  4 

13.277 

9.488 

p  =  12 

26.217 

21.026 
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Table  11.    Test  Results  of  Autocorrelation  and  Bivariate  Conditional  Normality 
in  the  BGARCH  Models  with  Positive  Definite  Parameterization 


Ii> 


Corn 


II 


2) 


Soybean 


1) 


111) 


Q(12) 


■St 

■ft 


10.84 
14.91 


10.84 
14.91 


20.23 
11.32 


20.23 
11.32 


Q2(12) 


^^t 

35.04 

35.21 

^\. 

52.41 

50.44 

Conditional 

Bivariate 

Normality 

Ti 

0.24 

0.33 

Tz 

-0.28 

-0.38 

Til 

4.35 

4.48 

T22 

3.54 

3.64 

c. 

34.44 

46.72 

39.60 
33.32 


40.59 
32.58 


-0.029 

.0000 

-0.13 

-0.23 

3.94 

3.80 

3.71 

3.60 

18.88 

15.76 

1) :  Refers  to  the  model  specification  without  drift  term, 

2):  Refers  to  the  model  specification  with  drift  term, 

T^:  Skewness,  i=l ,  2.  l=spot  returns,  2=futures  returns, 

T^i :  Kurtosis,  i=l ,  2.  l=spot  returns,  2=futures  returns, 

C3:  Bera-John's  LM  statistic,  x^(^) , 


Critical  values  of  x  (p)  ^^^e 


1% 

5% 

=     4 

13.277 

9.488 

=   12 

26.217 

21.026 
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Table  12.    Estimation  Results  of  the  Ordinary  Least  Square  Model,  Random 
Coefficient  Model,  and  Random  Coefficient  Autoregressive  Model, 
Corn 

OLS:     Rst   =  0.126  +  0.931  Rf^  +  e^  a^   -   39.413 

(0.357)   (0.029) 
r2    =   0.774,  D.W.  -  2.259, 
log  likelihood  =  -1012.61 

RCM:    Rst   =  0.358  +  ;3t  Rft  +  ^t 

(0.319) 
^t    =  0.987  +  /it 

(0.039) 
log  likelihood  =  -993.243 

RCMAR:   R^t   =  0.205  +  ^^  Rft  +  ^t 
(0.319) 
;9t    -   0.954  =  -0.758  (^^-i 

(0.029)   (0.123) 
log  likelihood  =  -985.393 

(  ) :  standard  error 


a2 

= 

20.040 
(2.530) 

^M 

0.1619 
(0.034) 

a^ 

= 

27.283 
(2.705) 

0.954)    +  /it 

^M 

= 

0.025 

(0.029) 

(0.015) 

Table  13.  Estimation  Results  of  the  Ordinary  Least  Square  Model,  Random 
Coefficient  Model,  and  Random  Coefficient  Autoregressive  Model, 
Soybean 

OLS:     Rgt   =  -0.128  +     0.897  Rft  +  Ct  a^   =   61.018 

(0.446)   (0.029) 
r2    =  0.703,  D.W.  =  1.634 
log  likelihood  =  -1080.58 

RCM:     Rst   =  -0.009  +  fi^   R^t  +  ^t 
(0.380) 
/9t   =  0.939  +  /it 
(0.046) 
log  likelihood  =  -1061.08 

RCMAR:   R^t   =  -0.0568+  p^   Rft  +  ^t 
(0.381) 
^t    -   0.934  =  0.280  (^t-i 
(0.050)   (0.134) 
log  likelihood  =  -1058.99 

(  ) :  standard  error 
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a2 

=  26.970 
(3.654) 

^M 

0.258 
(0.047) 

a^ 

=  30.070 
(3.860) 

0.934)    +  /it 

^M 

0.197 

(0.050) 

(0.044) 

Table  14.     In-Sample  Performance  of  Alternative  Hedging  Models  Compared 
to  No  Hedging 

Average  Percentage  Variance  Reduction 
from  the  No  Hedge  Position 


Corn 


•1) 


jlZ) 


Soybeans 


•1) 


II 


2) 


Constant  Hedging 
(Bivariate  Conditional) 


77.42 


-77.45 


-73.97 


-74.41 


Constant  Hedging(OLS) 


78.35 


78.78 


74.23 


74.40 


BGARCH 
(Constant  Corr) 


-78.75 


-78.21 


73.90 


-74.60 


BGARCH 
(Diagonal  Vec) 


79.42 


78.84 


-77.00 


76.99 


BGARCH 
(Positive  Def . ) 


78.92 


-78.94 


65.17 


-65.71 


RCMAR 


73.70 


73.22 


-73.23 


-73.55 


1) :    Refer  to  the  model  specification  without  drift  term 
2):   Refer  to  the  model  specification  with  drift  term 
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Table  15.    Out-of -Sample  Performance  of  Alternative  Hedging  Models  Compared 
to  No  Hedging 

Average  Percentage  Variance  Reduction 
from  the  No  Hedge  Position 


Corn 


•1) 


Il2) 


Soybeans 


1) 


II 


2) 


Constant  Hedging 
(Bivariate  Conditional) 


66.67 


52.22 


-80.84 


11.75 


Constant  Hedging(OLS) 


-67.70 


-52.23 


-80.80 


79.68 


BGARCH 
(Constant  Corr) 


■43.13 


-42.79 


-42.91 


17.34 


BGARCH 
(Diagonal  Vec) 


-69.61 


-54.30 


-85.69 


-84.91 


BGARCH 
(Positive  Def . ) 


-0.87 


6.49 


2.76 


2.77 


RCMAR 


-68.51 


-53.07 


17.07 


18.07 


1) :    Refer  to  the  model  specification  without  drift  term 
2):    Refer  to  the  model  specification  with  drift  term 
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